Mon. Not. R. Astron. Soc. OOP. [TlfTTI C2007') Printed 8 October 2009 (MN WT^ style file v2.2) 



Refined stellar, orbital and planetary parameters of the 
eccentric HAT-P-2 planetary system 



Andras Pal^'^'^*, Caspar A. Bakos^f, Cuillermo Torres^, Robert W. Noyes^, 
Debra A. Fischer^, John A. Johnson^, Cregory W. Henry^, R. Paul Butler^, 
Ceoffrey W. Marcy^, Andrew W. Howard^, Brigitta Sipocz^'^, 
David W. Latham^ and Cilbert A. Esquerdo^ 

Harvard-Smithsonian Center for Astrophysics, 60 Garden street, Cambridge, MA, 02138, USA 
^ Konkoly Observatory of the Hungarian Academy of Sciences, Konkoly Thege Miklos ut 15-17, Budapest, 1121, Hungary 
^ Department of Astronomy, Lordnd Eotvos University, Pdzmdny P. st. 1/A, Budapest, 1117, Hungary 
^ Department of Physics and Astronomy, San Francisco State University, San Francisco, CA, 94132, USA 
^ Institute for Astronomy, University of Hawaii, Honolulu, HI, 96822, USA 

^ Center of Excellence in Information Systems, Tennessee State University, Nashville, TN 37209, USA 

Department of Terrestrial Magnetism, Carnegie Institute of Washington, Washington DC, 20015, USA 
* Department of Astronomy, University of California, Berkeley, CA, 94720, USA 

Accepted . . . , received . . . ; in original form . . . 

ABSTRACT 

We present refined parameters for the extrasolar planetary system HAT-P-2 (also 
known as HD 147506), based on new radial velocity and photometric data. HAT- 
P-2b is a transiting extrasolar planet that exhibits an eccentric orbit. We present a 
detailed analysis of the planetary and stellar parameters, yielding consistent results 
for the mass and radius of the star, better constraints on the orbital eccentricity, 
and refined planetary parameters. The improved parameters for the host star are 
Mi, = 1.36 ± 0.04 M0 and i?* = 1.64 ± 0.08 i?©, while the planet has a mass of 
Mp = 9.09 ± 0.24 Mjup and radius of Rp ^ 1.16 ± 0.08i?jup. The refined transit 
epoch and period for the planet are E = 2, 454, 387.49375 ± 0.00074 (BJD) and P = 
5.6334729 ± 0.0000061 (days), and the orbital eccentricity and argument of periastron 
are e = 0.5171 ± 0.0033 and uj = 185.22° ± 0.95°. These orbital elements allow us to 
predict the timings of secondary eclipses with a reasonable accuracy of ~ 15 minutes. 
We also discuss the effects of this significant eccentricity including the characterization 
of the asymmetry in the transit light curve. Simple formulae are presented for the 
above, and these, in turn, can be used to constrain the orbital eccentricity using 
purely photometric data. These will be particularly useful for very high precision, 
space-borne observations of transiting planets. 

Key words: planetary systems — stars: fundamental parameters — stars: individual: 
HD 147506, HAT-P-2 - techniques: spectroscopic 



1 INTRODUCTION 

At the time of its discovery, HAT-P-2b was the longest pe- 
riod and most massive transiting extrasolar planet (TEP), 
and the only one known to exhibit an eccentric orbit 
IIBakos et al. 1 12007^ 1. In the following years, other TEPs 
have also been discovered with significant orbital eccen- 
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tricities and long periods: GJ 436b (iGillon et al. 2007 ) , 
HD 1 7156b (|Barbieri et al.ll2007l ). XO-3b (Ijohns-Krull et all 



2008h. and most notably HD 80606 (|Naef et all I2001I : 



Winn et al.ll2009bh . See http://exoplanet.eu for an up- 
to-date database for transiting extrasolar planets. 

The planetary companion to HAT-P-2 (HD 147506) 
was detected as a transiting obj ect during regular opera - 
tions of the HATNet telescopes l|Bakos et all I2OO2I . I2004D 
and the Wise HAT tel escope (WHAT, loca ted at the Wise 
Observatory, Israel; see lShporer et al"]|2006h . Approximately 
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26,000 individual photometric measurements of good signal- 
to-noise ratio (SNR) were gathered with the HATNet tele- 
scopes at the Fred Lawrence Whipple Observatory (FLWO, 
Arizona) and on Mauna Kea (Hawaii), and with the WHAT 
telescope. The planetary transit was followed up with the 
FLWO 1.2 m telescope and its KeplerCam detector. The 
planetary properties have been confirmed by radial velocity 
measurements and an analysis of the spectral line profiles. 
The lack of bisector span variations rules out the possibility 
that the photometric and spectroscopic signatures are due 
to a blended background eclipsing binary or a hierarchical 
system of three stars. 

The spin-orbit al ignment of the HAT -P-2( b) system was 
recent ly measured bv lWinn et al.l l|2007ah and lLoeillet et al] 
l|2008t ). Both studies reported an angle A between the pro- 
jections of the spin and orbital axes consistent with zero, 
within an uncertainty of ~ 10°. These results are particu- 
larly interesting because short period planets are thought to 
form at much larger distances and to then migrate inward. 
During this process, orbital ecce ntricity is tidally damped, 
yield ing an almost circular orbit l|D'Angelo. Lubow fc Batg 
I2OO6I ). Physical mechanisms such as Kozai interaction be- 
tween the transiting planet and an unknown massive com- 
panion on an inclined orbit could result in tight eccentric 
orbits (|Fabrvckv fc Tremaind l2007l : iTakeda. Kita fc Rasid 
I2OO8I ). However, in such a scenario, the spin-orbit align- 
ment as represented by A can be expected to be signifi- 
cantly larger. For instance, in the ca se of XO-3b, the re - 
ported alignments ar e A = 70° ±15° jHebrard et al.ll2008l ') 
and A = 37.3° ± 3.7° (|Winn et al.ll2009al . although there are 
indications of systematic observational effects). In multiple 
planetary systems, planet-planet scatte ring can also yield ec - 
centric and/or inclined orbits (see e.g. iFord fc Rasioll2008l ). 

The physical properties of the host star HAT-P-2 have 
been controversial, since different methods for stellar char- 
acterization have r esulted in stellar ra dii between ~ IARq 
and ~ 1.8 -R© (see iBakos et al. I l2007al ). Moreover, the true 
distance to the star has been uncertain in previous studies, 
with the Hipparcos-based distance being irreconcilable with 
the luminosity from stellar evolutionary models. 

In this paper we present new photometric and spec- 
troscopic observations of the planetary system HAT-P-2(b). 
The new photometric measurements significantly improve 
the light curve parameters, and therefore some of the stellar 
parameters are more accurately constrained. Our new radial 
velocity measurements yield significantly smaller uncertain- 
ties for the spectroscopic properties, including the orbital 
eccentricity, which have an impact also on the results of the 
stellar evolution modeling. In § [5] we summarize our photo- 
metric observations of this system, and in §[3] our new radial 
velocity measurements. The details of the analysis are dis- 
cussed in § 21 We summarize our results in § [S] 



2 PHOTOMETRIC OBSERVATIONS AND 
REDUCTIONS 

In the present analysis we make use of photometric data 
obtained with a variety of telescope/detector combinations, 
including the HATNet telescopes, the KeplerCam detector 
mounted on the FLWO 1.2 m telescope, the Nickel Im tele- 
scope at Lick Observatory on Mount Hamilton, California, 



and four of the automated photometric telescopes (APTs) 
at Fairborn Observatory in southern Arizona. The photo- 
me tric analysis of the HATNet data has been described 
by iBakos" et all l|2007ah . These HATNet data are shown in 
Fig. [1] with our new best-fit model superimposed (see § |4] 
for details on the light curve modeling). We observed the 
planetary transit on nine occasions: 2007 March 1^ (Sloan 
z band), 2007 April 21 (2), 2007 May 08 (z), 2007 June 22 
(z), 2008 March 24 (z), 2008 May 25 {z), 2008 July 26 {z), 
2009 April 28 (Stromgren h + y band), and 2009 May 15 
[b -\- y). These yielded 6 complete or nearly complete tran- 
sit light curves, and 3 partial events. One of these follow-up 
light curves (2007 April 21) was reported in the discovery 
paper. All of our individual high precision follow-up pho- 
tometry data are plotted in Fig. [51 along with our best-fit 
transit light curve model. The folded and binned light curve 
(computed only for the z-band observations) is displayed in 
Fig. El 

The frames taken with the KeplerCam detector were 
calibrated and reduced in the same way for the six nights 
at FLWO. For the calibrations we omitted saturated pixels, 
and applied standard procedures for bias, dark, and sky-flat 
corrections. 

Following the calibration, the detection of stars and 
the derivation of the astrometric solution was carried out 
in two steps. First, an initial astrometric transformation 
was derived using the ~ 50 brightest and non-saturated 
stars from each frame, and by using the 2MASS catalogue 
llSkrut skiel '20061) as a reference. We utilized the algorithm 
of jPa l & Bakog (|2006l ) with a second-order polynomial fit. 
The astrometric data from the 2MASS catalog were ob- 
tained from images with roughly the same SNR as ours. 
However, we expect significantly better precision from the 
FLWO 1.2 m owing to the larger number of individual ob- 
servations (by two orders of magnitude). Indeed, an internal 
catalog which was derived from the stellar centroids by reg- 
istering them to the same reference system has shown an 
internal precision of ~ 0.005 arc sec for the brighter stars, 
while the 2MASS catalog reports an uncertainty that is an 
order of magnitude larger: nearly ~ 0.06 arc sec. Therefore, 
in the second step of the astrometry, we used this new in- 
ternal catalog to derive the individual astrometric solutions 
for each frame, still using a second-order polynomial fit. We 
note here that this method also corrects for systematic er- 
rors in the photometry resulting from the proper motions of 
the stars, which have changed their position since the epoch 
of the 2MASS catalogue (~2000). 

Using the astrometric solutions above we performed 
aperture photometry on fixed centroids, employing a set of 
five apertures between 7.5 and 17.5 pixels in radius. The 
results of the aperture photometry were then transformed 
to the same instrumental magnitude system using a correc- 
tion to the spatial distortions and the differential extinction 
(the former depends on the celestial coordinates while the 
latter depends on the intrinsic colors of the stars). Both 
corrections were linear in the pixel coordinates and linear 
in the colors. Experience shows that significant correlations 
can occur between the instrumental magnitudes and some 



^ All of the dates are local (MST or HT) calendar dates for the 
first half of the night. 
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Fi gure 1. The folded HATNet light curve of HAT-P-2 (published 
in lBakoa et al.ir2007d) . showing the points only near the transit. 
The upper panel is superimposed with our best-fit model and the 
lower panel shows the residuals from the fit. See text for further 
details. 



of the external parameters of the light curves (such as the 
FWHM of the stars, and positions at the sub-pixel level). 
Ideally, one should detrend these correlations using only out- 
of-transit data (i.e., before ingress and after egress). Because 
of the lack of out-of-transit data, we instead carried out an 
external parameter decorrelation (EPD) simul taneous with 
the lig ht curve modeling (§ |4]) as described in iBakos et al.l 
(|2009l ) . After the simultaneous light curve modelling and de- 
trending, we chose the aperture for each night that yielded 
the smallest residual. In all cases this "best aperture" was 
neither the smallest nor the largest one from the set, con- 
firming the requirement to select a good aperture series. We 
note here that since all of the stars on the frames were well 
isolated, such choice of different radii for the apertures does 
not induce systematics related to variable blending of stars 
in different apertures. In addition, due to the high apparent 
brightness of HAT-P-2 and the comparison stars, the frames 
were acquired under a slightly extrafocal setting (in order to 
avoid saturation). This resulted in a different characteristic 
FWHM for each night. Thus, the optimal apertures yielding 
the highest SNR also have different radii for each night. Ad- 
ditional and more technical details abo ut the phot ometric 
reductions are discussed in Chapter 2 of EH (|2009b[ l. 

For the observations at Lick Observatory, we used the 
Nickel Direct Imaging Camera, which is a thinned Loral 
2048^ CCD with a 6.3' square field of view. We observed 
through a Gunn Z filter, and used 2x2 binning for an effec- 
tive pixel scale of 0'.'37 pixel" ^ . The exposure times were 25 s, 
with a readout and refresh time between exposures of 12 s. 
The conditions were clear for most of this transit with ~ I'.'O 
seeing. We defocused the images to draw out the exposure 
time while avoiding saturation for the target and reference 
stars. We applied the flat-field and bias calibrations, and 
determined the instrumental magnitude of HAT-P-2 using 
custom routines wr itte n in IDL as described previously by 
IWinn et al.] (l2007d ) and lJohnson et al.l (120081 ). We measured 
the flux of the target relative to two comparison stars using 
an aperture with a 17-pixel radius and a sky background 
annulus extending from 18 to 60 pixels. 

All four of the APTs at Fairborn Observatory have 



two channel photometers that measure t he Stromgre n 6 and 
Stromgren y count rates simultaneously (|Henrvll 19991 ). Since 
the Stromgren b and y bands are fairly close together and 
do not provide any useful color information for such shallow 
transits, we averaged the 6 and y differential magnitudes to 
create a {b + y) /2 "band pass" , which gives roughly a \/2 
improvement in precision. The comparison star for all of the 
APT observations is HD 145435. 



3 RADIAL VELOCITY OBSERVATIONS 



In the discovery paper for HAT-P-2b (|Bakos et al.ll2007al ') we 
reported thirteen individual radial velocity measurements 
from HIRES on the Keck I telescope, and ten radial veloc- 
ity measurements from the Hamilt on echelle spectrograph at 
the Lick Observatory l|Vogtlll9"87h . In the last year we have 
acquired 14 additional radial velocity measurements using 
the HIRES instrument on Keck. In the analysis we have 
incorporated as we ll the radial velocity data reported by 
iLoeillet et all (j2008h obtained with the OHP/SOPHIE spec- 
trograph. We use only their out-of-transit measurements, 
thereby avoiding the measurements affected by the Rossiter- 
McLaughlin effect. With these additional 8 observations, we 
have a total of 23 + 14 -|- 8 — 45 high-precision RV data 
points at hand for a refined analysis. 

In Table [T] we list all previously published RV measure- 
ments as well as our own new observations. These data are 
shown in Fig. |4l along with our best-fit model described be- 
low. 



4 ANALYSIS 

In this section we describe the analysis of the available pho- 
tometric and radial velocity data in order to determine the 
planetary parameters as accurately as possible. To model 
transit light curves taken in optical or near-infrared pho- 
tometric passbands, we include the effect of the stellar 
limb darkening. We h ave adopted the analytic formulae of 
iMandel fc Ago! (|2002D to model the flux decrease during 
transits under the assumption of a quadratic limb darken- 
ing law. Since the limb darkening coefficients are functions 
of the stellar atmospheric parameters (such as effective tem- 
perature Tcff, surface gravity logjr*, and metallicity) , the 
light curve analysis is preceded by an initial derivation of 
these parameters using the iodine-free template spectrum 
obtained with the HIRES instrument on Keck I. We em- 
plo yed the Spectrosc opy Made Easy software package (SME, 
see lValenti fc Piskun ov 1996), supported by the atomic line 



database of Valenti fc Fischerl l|2005l ). This analysis yields 
the Tcff, log g*, [Fe/H], and the projected rotational velocity 
V sin i. When all of these are free parameters, the initial SME 
analysis gives logg* = 4.22±0.14 (cgs), %f[ = 6290±110K, 
[Fe/H] = -f0.12±0.08, and ?;sini = 20.8 ± 0.2 kms^V The 
limb darkening coefficients were then derived for the z' , I, 
and (b + y) /2 phot ometric band s by interpo l ation, using the 
tables provided bv lClaretl l|2000l ) and lClaretl l|2004h . The ini- 
tial values for these coefficients were used in the subsequent 
global modeling of the data 14. ip . and also in refining the 
stellar parameters through a constraint on the mean stellar 
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Figure 2. Follow-up light curves of HAT-P-2. The top panel 
shows the z-band light curves acquired on 2007 March 18, 2007 
April 21, 2007 May 08, 2007 June 22, 2008 March 24, 2008 May 
25 and 2008 July 26; the respective transit sequence numbers are 
Nu- = -6, 0, +3, +11, +60, +71, and +82. The lower panel shows 
the Stromgren {b + y)/2 light curves, gathered on 2009 April 28 
and 2009 May 15, with transit sequence numbers A'^tr = +131 and 
+134. Our best-fit model is superimposed. See text for further 
details. 
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Figure 3. Folded and binned follow-up light curve of HAT-P- 
2, calculated from the seven individual 2-band events. The flux 
values at each point have been derived from ~35-50 individual 
measurements, and the bin size corresponds to a cadence of 3.6 
minutes (0.0025 days). The error bars are derived from the sta- 
tistical scatter of the points in each bin. Typical uncertainties are 
~ 0.4 mmag. 



density (see below). A second SME iteration was then per- 
formed with a fixed stellar surface gravity. The final limb 
darkening parameters are 7^^^' = 0.1419, 72^' = 0.3634, 



7f+^) = 0.4734, 7(''+^) 



= 0.2928, 



7^' 



= 0.1752, and 



0.3707. 



4.1 Light curve and radial velocity parameters 

The first step of the analysis is the determination of the 
light curve and radial velocity parameters. The parameters 
can be classified into three groups. The light curve parame- 
ters that are related to the physical properties of the plane- 
tary system are the transit epoch E, the period P, the frac- 
tional planetary radius p = Rp/Ri,, the impact parameter 
b, and the normalized semi-major axis a/Ri,. The physical 
radial velocity parameters are the RV semi- amplitude K, 
the orbital eccentricity e, and the argument of periastron 
uj. In the third group there are parameters that are not re- 
lated to the physical properties of the system, but are rather 
instrument-specific. These are the out-of-transit instrumen- 
tal magnitudes of the foUowup (and HATNet) light curves, 
and the zero-points 7Kcck, 7Lick, and 70HP of the three in- 
dividual RV data setfl. 

To minimize the correlation between the adjusted pa- 
rameters, we use a slightly different parameter set than that 
listed above. Instead of adjusting the epoch and period, we 
fitted the first and last available transit center times, T_i48 
and T+134. Here the indices denote the transit event number: 
the A'^tr = event was defined as the first complete follow-up 
light curve taken on 2007 April 21, the first available transit 



Since a synthetic stel lar spectrum was us ed as the reference 
in the reduction of the iLoeillet et al.l 1I2OO8I) data, 70HP is the 
actual barycentric radial velocity of the system. In the reductions 
of the Keck and Lick data we used an observed spectrum as the 
template, so the zero-points of these two sets are arbitrary and 
lack any real physical meaning. 
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Table 1. Complete list of relative radial velocity measurements 
for HAT-P-2 
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observation from the HATNet data was event A'tr = —148, 
and the last complete follow-up event (A'tr = +134) was ob- 
served on 2009 May 15. Note that if we assume the transit 
events are equally spaced in time, all of the transit centers 
available in the HATNet and follow-up photometry are con- 
strained by these two transit times. Similarly, instead of the 
eccentricity e and argument of periastron ui, we have used 
as adjustable parameters the Lagrangian orbital elements 
k = ecosu! and h = esintt;. These two quantities have the 
advantage of being uncorrelated for all practical purposes. 
Moreover, the radial velocity curve is an analytic function of 
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Figure 4. Radial velocity measurements for HAT-P-2 folded with 
the best-fit orbital period. Filled dots represent the OHP data, 
open circles show the the Lick/Hamilton, and the open boxes 
mark the Keck/HIRES observations. In the upper panel, all of 
these three RV data sets are shifted to zero mean barycentric 
velocity. The RV data are superimposed with our best-fit model. 
The lower panel shows the residuals from the best-fit. Note the 
different vertical scales in the two panels. The transit occurs at 
zero orbital phase. See text for further details. 

k and h even for cases where e — » (IPa il l2009al ). As is well 
known (jWinn et al.ll2007bl : [pill2008l ). the impact parame- 
ter b and a/7?* are also strongly correlated, especially for 
small v alues of p = R^/Ri,. Therefore, following the sugges- 
tion bv lBakos et al.l lj2007bl ). we have chosen the parameters 
C,/R* and for fitting instead of a/R^ and h, where C/^* 
is related to a/_R* as 

a \ 27r 

'r 



r^ 



p VT 



(1) 



-62 l + h 

is related to the transit duration as 



The quantity (^/R 
Tdur = 2(c^/i?*)~^, the duration here being defined between 
the time instants when the center of the planet crosses the 
limb of the star inwards and outwards, respectively. 

The actual flux decrease caused by the transiting planet 
can be estimated from the projected radial distance between 
the center of the planet and the center of the star d (nor- 
malized to -R*). Fo r circular orbits the t ime dependence of 
d is trivial (see e.g. lMandel fc Ago]||200j ). For eccentric or- 
bits, it is necessary to use a precise paramctrization of d as 
a function of time. As was shown by Pal (2008_), d can be 
expressed in a second order approximation as 



d^ = (1 



^) {Atf + b\ 



(2) 



where At is the time between the actual transit time and 
the RV-based transit center. Here the RV-based transit cen- 
ter is defined when the planet reaches its maximal tangen- 
tial velocity during the transit. Throughout this paper we 
give the ephemeris for the RV-based transit centers and de- 



© 2007 RAS, MNRAS OOO.mfTTI 



6 A. Pal et al. 



note these simply by Tc. Although the tangential velocity 
cannot be measured directly, the RV-based transit center is 
constrained purely by the radial velocity data, without re- 
quiring any prior knowledge of the transit geometrjjf]. For 
eccentric orbits the impact parameter 6 is related to the or- 
bital inclination i by 



h 



(3) 



In order to have a better description of the transit light 
curve, we used a higher order expansion in the d(At) 
function (Eq. [2]). For circular orbits, such an expansion is 
straightforward. To derive the expansion for elliptical or- 
bits, we employed the method of Lie-integration which gives 
the solution of any ordinary differential equation (here, the 
equations for the two-body problem) in a recursive series for 
the Taylor expansion with respect to the independent vari- 
able (here, the time). By substituting the initial conditions 
for a body of which spatial coordinates are written as func- 
tions of th e orbit al elements, using equations (C1)-(C8) of 
iPal fc Siilil \2m% one can derive that the normalized pro- 
jected distance d up to fourth order is: 



d' = 
where 

Q = 

and 
R = 



2 



+ 



l + h 



p2-)3/: 



-k. 



(4) 



(5) 



(6) 



Here n = 2-k / P is the mean motion, and <^ is defined as 
ip — nAt. For circular orbits, Q — 1 and R — 0, and for 
small eccentricities (e -C 1), Q ~ 1 -I- 3/i and R^ k. 



4.2 Joint fit 

Given the physical model parametrized above, we performed 
a simultaneous fit of all of the light curve and radial velocity 
data. We used Eq. Q to model the light curves, where the 
parameters Q and R were derived from the actual values of 
k and h, using equations Eq. ^ and Eq. To find the 
best-fit values for the p arameters we emp loyed the downhill 
simplex algorithm (see iPress et al.l Il992l ) and we used the 
method of refitting to synthetic data sets to infer the proba- 
bility distribution for the adjusted values. In order to char- 
acterize the effects of red noise properly, the mock datasets 
in this bootstrap method were generated by perturbing ran- 
domly only the phases in the Fourier spectrum of the residu- 
als. The final results of the fit were T-ms = 2453379. 10210± 
0.00121, r+134 = 2454967.74146 ± 0.00093, K = 983.9 ± 
17.2 ms-\ k = -0.5152 ± 0.0036, h = -0.0441 ± 0.0084, 



Rp/R^, = p = 0.07227 ± 0.00061, = 0.156 ± 0.074, 
C/i?* = 12.147 ± 0.046 day-S 7Kcck = 316.0 ± 6.0ms-\ 
7Lick = 88.9±10.4ms~\ and 70HP = -19860.5±10.2 ms"\ 
The uncertainties of the out-of-transit magnitudes were in 
the range 6-21 x 10~^ mag for the follow-up light curves, and 
16 X 10"^ mag for the HATNet dat<|3. The fit resuhed in a 
reduced value of 0.992. As described in the following sub- 
section, the resulting distributions of parameters have been 
used subsequently as inputs for the stellar evolution model- 
ing. 



4.3 Effects of the orbital eccentricity on the 
transit 

In this section we summarize how the orbital eccentricity af- 
fects the shape of the transit light curve. The leading-order 
correction term in Eq. Q in —2b^R(p, is related to the 
time lag betw een the phot ometric and RV-based transit cen- 
ters (see also lKopailll959l ). The photometric transit center, 
denoted Tc,phot , is defined halfway between the instants 
when the center of the planet crosses the limb of the star 
inward and outward. It is easy to show by solving the equa- 
tion d{if) = 1, yielding two solutions {ipi and V'e), that this 
phase lag is: 



b^R 



a 



62) - (Q - i?2)b2 

b^k 



{i + h)Vi 



(7) 
(8) 

(9) 



which can result in a time lag of several minutes. For in- 
stance, in the case of HAT-P-2b, Tc.phot — 7c,rv ~ n~ A(p — 
1.6 ± 0.9 minutes. 

In Eq. the third order terms in if describe the asym- 
metry between the slopes of the ingress and egress parts of 
the light curve. For other aspect s of li ght curve asymme- 
tries, see iLoebl l|2005l ) and lBarned (|2007l ). In cases where no 
constraints on the orbital eccentricity are available (such as 
when there are no RV measurements), one cannot treat the 
parameters R and Q as independent since the photometric 
transit center and R have an exceptionally high correlation. 
However, if we assume a simpler model function, with only 
third order terms in ip with fitted coefficients present, i.e. 



d' = b' 



1 2 1^3 

1-ifi - -Cip 



3^ 



(10) 



these will yield a non-zero value for the C coefficient for 
asymmetric light curves. In the case of HAT-P-2b, the de- 
rived values for Q and 7? are Q = 2.204 ± 0.074 and 
R = —0.784 ± 0.015 (obtained from the values of k and 
h; see § 14. 2p . Therefore, the coefficient for the third-order 
term in ifi will be QR = -1.73±0.09. Using Eq. CO}, for an 
"ideal" light curve (with similar parameters of k, h, (,/Ri, 



^ In other words, predictions can only be made for the RV-based 
transit center in the cases where the planet was discovered by a 
radial velocity survey and initially there are no further constraints 
on the geometry of the system, notably its impact parameter. 



* Note that these small uncertainties reflect only the uncertain- 
ties of the instrumental magnitudes, and not the intrinsic magni- 
tudes in some absolute photometric system. 
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Table 2. Stellar parameters for HAT-P-2. 



Parameter 


Value 


Source 


TeS (K) 


6290 ± 60 


O TV /fT~^P\ 

SME 


[Fe/H] 


+0.14 ± 0.08 


SME 


logs* (cgs) 


4.16 ±0.03 


SME 


t!sini (kms"-*^) 


20.8 ± 0.3 


SME 


(M©) 


1.36 ±0.04 


Y^+LC+SME'^ 


R* (Rq) 




Y2+LC+SME 


logs* (cgs) 


4.138 ±0.035 


Y^+LC+SME 


L* (Lo) 


q 70+0.48 
■J- '"-0.38 


Y^+LC+SME 


My (mag) 


3.31 ±0.13 


Y^+LC+SME 


Age (Gyr) 


2.6 ±0.5 


Y^+LC+SME 


Distance (pc) 


119 ±8 


Y^+LC+SME 



SME = "Spectrosco py Made Easy" package for analysis of high- 
resolution spectra bv lValenti fe Piskunovl lll996l). See text . 
^ Y^+LC+SME = Yonsei-Yale isochrones liYi et al.ll200lh . light 
curve parameters, and SME results. 



and as for HAT-P-2b), the best fit value for C will be 
C = —2.23, which is close to the value of QR ~ —1.73. The 
difference between the best fit value of C and the fiducial 
value of QR is explained by the fact that in Eq. (|10p we 
adjusted the coefficient for the third order term in (p that 
causes the asymmetry in the light curve, and therefore the 
corrections in the lower-order terms (such as —2R, Q — R^, 
and Q/3 in Eq. |4ll have been neglected. 

Although this asymmetry can in principle be measured 
directly (without leading to any degeneracy between the fit 
parameters) , in practice one needs extreme photometric pre- 
cision to obtain a significant detection for a non-zero C pa- 
rameter. Assuming a photometric time series for a single 
transit of HAT-P-2b with 5 sec cadence where each individ- 
ual measurement has a photometric error of 0.01mmag(!), 
the uncertainty in C will be ±0.47, equivalent to a 5-a detec- 
tion of the light curve asymmetry. This detection would be 
difficult with ground-based instrumentation. For example, 
for a 1-a detection one would need to achieve a photometric 
precision of 0.05 mmag at the same cadence, as suming purely 
white noise). Space missions such as Kepler (|Borucki et al.l 
I2OO7I ) will be able to detect orbital eccentricity of other plan- 
ets relying only on transit photometry. 

4.4 Stellar parameters 

As pointed out by ISozzetti et all (|2007t ). the ratio a/J?* is 
a more effective luminosity indicator than the spectroscop- 
ically determined stellar surface gravity. In the cases where 
the mass of the transiting planet is negligible, the mean stel- 
lar density is 

The normalized semi-major axis a/Ri, can be obtained from 
the transit light curve model parameters, the orbital eccen- 
tricity, and the argument of periastron (see Eq. H}). 

Since HAT-P-2b is quite a massive planet, (Mp/M* ~ 
0.01), relation (|lip requires a small but significant cor- 
rection, which a lso depends on observable quantities (see 
iPal et al.|[2008bl , for more details). For HAT-P-2b this cor- 
rection is not negligible because Mp/M* is comparable to 



the typical rel ative uncertainties in the light curve parame- 
ters. Following |Pal et al.l l|2008al ) the density of the star can 
be written as 



So 

where both po and Eq are observables, namely, 
3 

Po = 



GP2 \ R. 



So = 



2PG sin i 



(12) 

(13) 
(14) 



In Eq. p2|l the only unknown quantity is the radius of the 
star, which can be derived using a stellar evolution model, 
and it depends on a luminosity indicatoi[f|, the effective tem- 
perature TcfT (obtained from the SME analysis), and the 
chemical composition [Fe/H]. Therefore, one can write 



7?. = 7?.(p.,Teff,[Fe/H]). 



(15) 



Since both T^s and [Fe/H] are known, we may solve for the 
two unknowns in Eq. (|12|l and Eq. (|15|l . Note that in order 
to solve Eq. (|15p . supposing its parameters are known in ad- 
vance, one needs to make use of a certain stellar evolution 
model. Such models are available only in tabulated form, 
and therefore the solution of the equation requires the in- 
version of the interpolating function on the tabulated data. 
Thus, Eq. (|15|l is only a symbolical notation for the algo- 
rithm which provides the solution. Moreover, if the star is 
evolved, the isochrones and/or evolutionary tracks for the 
stellar models can intersect each other, resulting in an am- 
biguous solution (i.e., one no longer has a "function" , strictly 
speaking). For HAT-P-2, however, the solution of Eq. (fT5)l 
is definite since the host star is a relatively unevolved main 
sequence star. To obtain the physical parameters (such as 
the st ellar radius) we used the evolution models of lYi et ahl 
(|200ll ). and interpolated the valu es of p<-, Teff and [Fe/H ] 
using the interpolator provided bv lDemargue et al.1 (120041 ) ■ 
The procedure described above has been applied to all 
of the parameters in t he input set in a complete Monte- 
Carlo fashion (see also |Pal et al.|[2008ah . where the values 
of Po have been derived from the values of a/i?* and the 
orbital period P using Eq. (|13|) . while the values for Toff 
and [Fe/H] have been drawn from Gaussian distributions 
with the mean and standard deviation of the first SME re- 
sults (Teff = 6290 ± 110 K and [Fe/H] = ±0.12 ± 0.08). 
This step produced the probability distribution of the phys- 
ical stellar parameters, including the surface gravity. The 
value and associated uncertainty for that particular quan- 
tity is logg* = 4.16 ± 0.04 (cgs), which is slightly smaller 
than the result from the SME analysis. To avoid systematic 
errors in Toff and [Fe/H] stemming from their correlation 
with the spectroscopically determined (and usually weakly 
constrained) logg*, we repeated the SME analysis by fix- 
ing the value of logy* to the above value from the mod- 
eling. This second SME run gave T^s = 6290 ± 60 K and 
[Fe/H] = ±0.14 ± 0.08. We then updated the values for the 
limb darkening parameters, and repeated the simultaneous 
light curve and radial velocity fit. The results of this fit were 

^ In practice this is either the surface gravity, the density of the 
star, or the absolute magnitude (if a parallax is available). 
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then used to repeat the stellar evolution modeling, which 
yielded among other parameters logg* = 4.138±0.035 (cgs). 
The change compared to the previous iteration is small 
enough that no further iterations were necessary. Our use 
here of this classic treatment of error propagation instead 
of a Bayesian approach in order to derive the final stellar 
parameters is essentially determined by the functionalities 
of the SME package. In view of the fact that the surface 
gravity from the stellar evolution modeling (constrained by 
the photometric and RV data) has a significantly smaller 
uncertainty than the value delivered by the SME analysis, 
we believe this kind of iterative solution and the method of 
error estimation are adequate. The stellar parameters are 
summarized in Table (2] and the light curve and radial ve- 
locity parameters are listed in the top two blocks of Table |31 

4.5 Planetary parameters 

In the two previous steps of the analysis we determined the 
light curve, radial velocity curve, and stellar parameters. In 
order to obtain the planetary parameters, we combined the 
two Monte-Carlo data sets providing probability distribu- 
tions for all quantities in a consistent way. For example, the 
mass of the planet is calculated using 



2-K KVT 



Gsini 



(16) 



where the values for the period P, RV semi-amplitude K, 
eccentricity e, inclination i, and normalized semi-major axis 
a/Rk were taken from the results of the light curve and 
RV fit, while the values for i?* were taken from the cor- 
responding stellar parameter distribution. From the distri- 
bution of the planetary parameters, we obtained the mean 
values and uncertainties. We derived Mp — 9.09 ± 0.24 Mjup 



for the planetary mass and Rp = 1.157 



+0.073 



R]up for the 



radius, with a correlation coefficient of C{Mp, Rp) = 0.68 
between these two parameters. The planetary parameters 
are summarized in the thir d block of Table [3] Compared 
with the values reported bv iBakos" et al.1 l|2007ah . the mass 
of the planet has not changed significantly (from Alp = 
9.04 ± 0.50 Mjup), but the uncertainty is now smaller by 
a factor of two. The new estimate of the planetary radius 
is larger by roughly 2-a, while its uncertainty is similar 
or shghtly smaller than before {Rp = 0.982t|5 ?i55 Rjup for 
iBakos et al]|2007ah . 

The surface temperature of the planet is poorly con- 
strained because of the lack of knowledge about redistribu- 
tion of the incoming stellar flux or the effects of significant 
orbital eccentricity. Assuming complete heat redistribution, 
the surface temperature can be estimated by time averaging 
the incoming flux, which varies as 1/r^ — {1 — e cos E)~^ 
due to the orbital eccentricity. The time average of 1/r^ is 



dt 



1 

2^ 



dM 

'■2(M)' 



(17) 



where M is the mean anomaly of the planet. Since r = a(l — 
e cos E) and dM = (1 — e cos E)dE, where E is the eccentric 
anomaly, the above integral can be calculated analytically 
and the result is 



1 



(18) 



Table 3. Spectroscopic and light curve solutions for HAT-P-2, 
and inferred planetary parameters. 



Parameter 



Value 



Light curve parameters 
P (days) 

E (BJD - 2,400,000) 

Ti4 (days)^ 

Ti2 = T34 (days)^ 

a/R^, 

Rp/R-k 

b = (a/i?*) cosi(l - e2)/(l -l-h) 
i (deg) 

Spectroscopic (RV) parameters 
K (ms~i) 
k = ecosw 
h = esintj 



Planetary parameters 
Mp (Mjup) 

Rp (Rjup) 
C{Mp,Rp) 
Pp (g cm-3) 
a (AU) 
logffp (cgs) 
Teff (K) 
Secondary eclipse 



5.6334729 ± 0.0000061 
54, 387.49375 ± 0.00074 
0.1787 ±0.0013 



0141+°-°"i'5 

□+0.39 
-0.41 



8.99: 
0.07227 ± 0.00061 



0.395in^g 
hi. 12 



86°. 72"' 



0.87 



Esc 
^14, sec 



(BJD - 2,400,000) 
(days) 



983.9 ± 17.2 
-0.5152 ±0.0036 
-0.0441 ± 0.0084 
0.5171 ±0.0033 
185.22° ±0.95° 

9.09 ± 0.24 

1 1C7+0.073 

1.10(_Q Qg2 

0.68 
7.29 ± 1.12 
0.06878 ± 0.00068 

4.226 ± 0.043 
1540 ± 30 (see ^) 

0.1868 ± 0.0019 
54,388.546 ± 0.011 
0.1650 ±0.0034 



^ T14: total transit duration, time between first and last contact; 
T12 = T34: ingress/egress time, time between first and second, or 
third and fourth contact. 

^ This effective temperature assumes uniform heat redistribution, 
while the irradiance is averaged over the entire orbit. See text for 
further details about the issue of the planetary surface tempera- 
ture. 

Using this time-averaged weight for the incoming flux, we 
derived Tp — 1540 ± 30 K. However, the planet surface tem- 
perature would be ~ 2975 K on the dayside during perias- 
tron assuming no heat redistribution, while the equilibrium 
temperature would be only ~ 1190K at apastron. Thus, we 
conclude that the surface temperature can vary by a factor 
of ~ 3, depending on the actual atmospheric dynamics. 



4.6 Photometric parameters and the distance of 
the system 

The measured col or index of the sta r as reported in the 
TASS catalogue (|Droege et al.l l2006t ) is (V - 7)tass = 
0.55±0.06, which is in excellent agreement with the result of 
{V — I)yy ~ 0.552 ± 0.016 we obtain from the stellar evolu- 
tion modeling (see § I4.4|l . The models also provide the abso- 
lute visual magnitude of the star as My ~ 3.31±0.13, which 
gives a distance modulus of Vtass — My ~ 5.39 ±0.13 corre- 
sponding to a distance of 119±8 pc, assuming no interstellar 
extinction. This distance estimate is intermediate between 
the values inferred from the trigonometric parallax in the 
original Hipparcos catalog (tthip = 7.39 ± 0.88 mas, corre - 
sponding to a distance of 135± 18 pc: |Perrvman et al.lll997l ). 
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Figure 5. Observational constraints for HAT-P-2 compared with 
stellar evolution calculations from the Yonsei-Yale models, repre- 
sented by isochrones for [Fe/H] = +0.14 between 0.5 and 5.5 Gyr, 
in steps of 0.5 Gyr. The luminosities on the vertical axis are ren- 
dered in two ways: as absolute visual magnitudes My in the top 
panel, and with the ratio a/Ri, as a proxy in the lower panel. The 
effective temperature along with the absolute magnitudes inferred 
from the apparent brightness in the TASS catalog and the origi- 
nal and revised Hipparcos parallaxes are shown in the top panel 
with the corresponding l-cr and 2-cr confidence ellipsoids (upper 
ellipsoid for the original Hipparcos reductions, lower for the re- 
vision). The diamond represents the value of My derived from 
our best- fit stellar evolution models using a/Ri, as a constraing 
on the luminosity. In the lower panel we show the confidence el- 
lipsoids for the temperature and our estimate of a/Ri, from the 
light curve. 



and in the revised r eduction of the original Hipparcos obser- 
vations by Ivan Leeu won (2007a, b) (tthip = 10. 14± 0.73 mas, 
equivalent to a distance of 99 ± 7pc). In Fig. [5] the model 
isochrones are shown for the measured metallicity of HAT- 
P-2 against the measured effective temperature and two sets 
of luminosity constraints: those provided by the estimates of 



44 




58 f ^ 

0.0 0.1 0.2 0.3 0.4 0.5 

Mean motion (d"') of a possible second planet 



Figure 6. Unbiased of ^ three-body Keplerian + circular fit 
to the RV observations, where the mean motion of a possible 
secondary companion has been varied between and 0.5 d~^. 

the Hipparcos distance (original, and revised) together with 
the TASS apparent magnitudes (top panel), and the con- 
straint from the stellar density inferred from the light curve 
(bottom panel). We note in passing that the distance derived 
using the near-infrared 2MASS photometry agrees well with 
the distance that relies on the TASS optical magnitudes. 

4.7 Limits on the presence of a second companion 

In this section we discuss limits on the presence of an ad- 
ditional planet in this system. We performed two types of 
tests. In both of these tests we have fitted the RV semi- 
amplitude K, the Lagrangian orbital elements {k,h), the 
three velocity zero-points (7Keck, 7Lick and 7ohp), and the 
additional terms required by the respective test methods 
(drift coefficients or orbital amplitudes). In these fits, the 
orbital epoch E and period P of HAT-P-2b have been kept 
fixed at the values yielded by the joint photometric and RV 
fit. This is a plausible assumption since without the con- 
straints given by the photometry, the best fit epoch and 
period would be Erv = 2454342.455 ± 0.016 (BJD) and 
P — 5. 6337± 0.0016, i.e., the uncertainties would be roughly 
20-25 times larger. 

In the first test, a linear, quadratic and cubic polyno- 
mial were added to the radial velocity model functions in 
addition to the 7 zero-point velocities. Fitting a linear trend 
yielded a drift of Guncar = —21.2 ± 12.1 ms~^ yr~^, with 
the decreasing from 52.1 to 48.5 (note that in this test 
the effective number of degrees of freedom is 45 — 7 = 38) . 
Therefore, both the decrease in the residuals and the rela- 
tive uncertainty of Gunoar suggest a noticeable but not very 
significant linear drift on the timescale of the observations 
(that is, approximately, 1.7years). The additional quadratic 
and cubic terms do not yield a significant decrease in the 
unbiased residuals. 

In the second test, we extended the system configura- 
tion with an additional planet orbiting the star on a circular 
orbit. The orbital phase and the semi-amplitude of this addi- 
tional companion were fitted simultaneously with the Keple- 
rian orbital elements of HAT-P-2b, while the mean motion 
of the second companion was varied between 712 — 0.001 
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and n2 — 0.5 d ~ 0.45nHAT-p-2 with a step size of 
An = 0.001 d"^ < (1.7yr)~\ As can be seen in Fig. |6l 
no significant detection of a possible secondary companion 
can be confirmed. 



4.8 Secondary eclipse timings 

The improved orbital eccentricity and argument of perias- 
tron allow us to estimate the time of the possible occulta- 
tions. For small orbital eccentricities, the offset of the sec- 
ondary eclipse from phase 0.5 is proportional to k = ecosui 
(see, e.g., Charbonncau ot al. 2005,). However, in the case of 
larger eccentricities as in HAT-P-2b, this linear approxima- 
tion can no longer be applied. The appropriate formula for 
arbitrary eccentricities can be calculated as the difference 
between the mean orbital longitudes at secondary eclipse 
(Ascc) and at transit (Apri), that is, 

2kJ 



than the measured radius of l.lGlg Qg 



A, 



= 7r + 



+ arg 



/l2 



7 2 2 

k e 



(19) 



2k^ 



(1 + J)2 1 + J 



2k - 



2e^k 



l + J 



where J — \/l — e^. It is easy to see that the expansion of 
Eq. lUll yields 

Ascc - Apri ^ TT + 4fc (20) 

fo r \k\ <C 1 and \h\ <^ 1, an d this is equivalent to Eq. (3) 
of lCharbonneau et al.l (|2005l '). In the case of HAT-P-2b, we 
find that secondary eclipses occur at the orbital phase of 
flisec = (Ascc-Apri)/(27r) = 0.1868±0.0019, i.e., Iday 1 hour 
and 17 minutes (± 15 minutes) after the transit events. 



5 DISCUSSION 

In this work we have presented refined planetary, stellar and 
orbital parameters for the HAT-P-2(b) transiting extrasolar 
planetary system based on a full modeling of new and exist- 
ing data. These data consist of previously published radial- 
velocity measurements along with new spectroscopic obser- 
vations, and a new set of high-precision photometric obser- 
vations of a number of transit events. The refined parameters 
have uncertainties smaller by a factor of ~ 2 in the planetary 
parameters and a factor of ~3-4 in t he orbital pa r ameter s 
than the previously reported values of iBakos et al.l (l2007al ). 
We note that the density of the planet as determined here, 
Pp = 7.29 ± 1.12 gcm~'^, is significantly smal ler than the 
value P p,b2007 ~ 11.9ti'.6 gcm"^ inferred bv iBakos et al.l 
(|2007al ). and the new uncertainty is significantly smaller as 
well. Our analysis does not rely on the distance of the sys- 
tem, i.e., we have not made use of the absolute magnitude as 
a luminosity indicator. Instead, our stellar evolution model- 
ing is based on the density of the star, which is a proxy for 
luminosity and can be determined to high precision directly 
from photometric and RV observations. A comparison of the 
distance of the system as derived from the model absolute 
magnitude with the Hipparcos determination (original, and 
revised) shows that our (density-based) estimate is interme- 
diate between the two astrometric determinations. 

The ze ro ins olation planetary isochrones of 
iBaraffe et al.l (|2003l ) give an expected radius of 
i?p,Baraffc03 = 1-02 ± 0.02_Rjup, which is slightly smaller 



-„uu^Jup- The work of 

iFortnev. Marlev fc Barne3 |20o3) takes into account not 
only the evolutionary age and the total mass of the planet, 
but the incident stellar fiux and the mass of the planet's core 
as well. By scaling the semi-major axis of HAT-P-2b to one 
that yields the same incident fiux from a solar-type star on 
a circular orbit, taking into account both the luminosity of 
the star and the correction for the orbital eccentricity given 
by Eq. (0, we derived a' = 0.033 ± 0.003 AU. Using this 
scaled semi major axis, the interpolation based on the tables 
provided by iFortnev. Marlev Sz Barne^ l|2007h yields radii 
between -Rp,Fortney,o = 1.142 ±0.003 Tijup (core-less planets) 
and -Rp.Fortncy.ioo ~ 1 . 1 1 1 ± 0.003 -Rjup (core-dominated 
planets, with a core of Afp^corc ~ 100 Mq). Although these 
values agree nicely with our value of Rp — 1.1571q 

R Jup , 

the relatively laxge uncertainty of Rp precludes any further 
conclusions about the size of the plan e t's co re. Recent 
models by iBaraffe. Chabrier fc BarmanI (|2008l ) also give 
the radius of the planet as a function of evolutionary age, 
metal enrichment, and an optional insolation equivalent to 
a' = 0.045 AU. Using this latter insolation, their models 
yield iip,Baraffe08,o.02 = 1 .055 ± 0.006 i?jup (for metal poor, 
Z = 0.02 planets) and i?p,Baraffe08,o.io = 1.008 ± 0.006 i?jup 
(for more metal rich, Z = 0.10 planets). These values 
are slightly smaller than the actual radius of HAT- 
P-2b. However, the actual insolation of HAT-P-2b is 
roughly two times larger than the insolation implied 
by a' = 0.045 AU. Since t he p lanetary radius from 
IBaraffe. Chabrier fc BarmanI (120081 ) for zero insola- 
tion gives ^?^°L,affc08,o.02 = 1-009 ± 0.006 Ejup and 
^p°Baraffc08 10 ~ 0.975 ± 0.006 i?jup for metal enrichments 
of Z = 0.02 and Z — 0.10, respectively, an extrapolation 
for a two times larger insolation would put the expected 
planetary radius in the r ange of ^ 1.10 -Rj up. This is consis- 
tent with the models of IFortnev. Marlev fc Barne j (|2007l ) 
as well as with the mea surements. However, as discussed 
earlier in the case of the IFortnev. Marlev fc Barne^ (|2007| ) 
models, the uncertainty in Rp does not allow us to properly 
constrain the metal enrichment for the recent Baraffe 
models. 
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